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SUMMARY 


Flight  tests  with  the  NASA  Terminal  Configured  Vehicle  (TCV)  Boeing  737 
were  made  to  investigate  the  use  of  back-azimuth  signals  to  estimate  airplane 
position.  Position  estimates  were  established  by  three  update  methods:  a 
combination  of  microwave  landing  system  (MLS)  back-azimuth  and  distance¬ 
measuring-equipment  (DME)  signals,  MLS  back-azimuth  signals  alone,  and  dual  DME 
signals.  Position-estimate  error  was  obtained  by  comparison  of  the  position 
estimate  from  the  airplane  navigation  system  with  the  tracking  information  from 
a  ground-based  radar  system. 

The  results  of  these  tests  showed  that  the  most  accurate  position  esti¬ 
mates  were  obtained  with  a  combination  of  back -azimuth  and  DME  signals.  The 
next  most  accurate  estimates  were  made  with  only  back -azimuth  signals.  Of  the 
three  update  modes  tested,  the  least  accurate  estimates  were  made  with  dual 
DME  signals. 

Analysis  of  the  position-estimate  error  data  showed  that  the  component  of 
error  due  to  the  DME  signal  input  was  larger  than  the  component  due  to  the  MLS 
back -azimuth  signal  input.  The  use  of  a  precision  DME  signal  input  would  be 
expected  to  reduce  position-estimate  error. 


INTRODUCTION 

The  NASA  Terminal  Configured  Vehicle  (TCV)  previously  demonstrated  the 
ability  to  navigate  along  a  curved  path  approach  to  landing  (ref.  1)  with  an 
onboard  navigation  system  utilizing  position  information  from  the  time- 
referenced  scanning-beam  (TRSB)  microwave  landing  system  (MLS) .  During 
operations  in  the  MLS  environment,  azimuth  and  elevation-angle  information  was 
obtained  from  two  scanning  beams,  and  distance  information  was  provided  by  MLS 
precision  distance-measuring  equipment  (PDME)  for  use  in  the  onboard  navigation 
and  automatic  landing  systems.  The  TRSB  MLS  provides  azimuth  information 
within  60°  of  each  side  of  the  runway  center  line,  elevation-angle  information 
up  to  20°,  and  PDME  range  information  up  to  37  040  m  (20  n.  mi.)  from  the  MLS 
antennas . 

Missed-approach  navigation  capability  was  available  from  the  MLS  system 
in  the  form  of  back-azimuth  and  PDME  information.  Hence,  to  fully  utilize  the 
TRSB  MLS,  algorithms  using  back-azimuth  information  were  incorporated  into  the 
TCV  airplane  navigation  computer  for  use  in  generating  missed-approach  guidance. 
Flight  tests  then  were  conducted  to  determine  the  navigation  accuracy  from  MLS 
back-azimuth  signals. 

The  purpose  of  this  report  is  to  describe  the  equations  and  logic  used  to 
generate  a  navigation  position  estimate  in  the  MLS  back-azimuth  signal  environ¬ 
ment  and  to  document  the  error  in  the  navigation  position  estimate.  The  equa¬ 
tions  described  are  used  to  calculate  position  difference  components  from 
(1)  MLS  back-azimuth  bearing  and  range  information  from  arbitrarily  located 


DME  stations,  and  (2)  MLS  back-azimuth  bearing  information  alone.  A  summary 
of  the  TCV  position-estimate  update  process  is  also  described.  The  navigation 
position-estimate  error  calculated  from  flight  data  and  radar  tracking  infor¬ 
mation  is  analyzed  to  determine  the  relative  errors  due  to  DME  and  back-azimuth 
signal  inputs.  The  position-estimate  error  data  using  the  MLS  inputs  are  also 
compared  with  error  data  obtained  during  dual  DME  updates,  the  primary  position- 
estimate  mode. 

SYMBOLS 

A  distance  between  the  DME  and  MLS  back-azimuth  antenna,  n.  mi. 

a  east  component  of  A,  n.  mi. 

b  north  component  of  A,  n.  mi. 

D  distance  from  DME  to  airplane  corrected  for  slant  range,  n.  mi. 

D'  slant  range  distance  from  DME  to  airplane,  n.  mi. 

DP  position  difference  vector 

DPp  component  of  position  difference  vector  DP  perpendicular  to  the 

runway  center  line 

DPp  magnitude  of  vector  DPp,  n.  mi. 

ep  magnitude  of  the  component  of  position-estimate  error  perpendicular 

to  the  back-azimuth  bearing  on  which  the  airplane  is  located,  m 

eR  magnitude  of  the  component  of  position-estimate  error  parallel  to 

the  back-azimuth  bearing  on  which  the  airplane  is  located,  m 

eg  azimuth  bearing-error  angle,  deg 

F  ellipticity  constant,  0.003367 

ha/p  altitude  of  airplane  above  mean  sea  level,  m 

hoME  altitude  of  DME  antenna  above  mean  sea  level,  m 

h0  altitude  of  MLS  back-azimuth  antenna  above  mean  sea  level,  m 

i,j  unit  coordinate  vectors 

Ki  position  update  gain 

K2  velocity  update  gain,  sec”1 

L  outer  radial  limit  of  back-azimuth  volumetric  check,  m  (n.  mi.) 
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M  inner  radial  limit  of  back-azimuth  volumetric  check,  m 

N,E  axes  of  orthogonal  coordinate  system  oriented  toward  true  north 

P  angle  formed  by  the  vector  Z,n  and  a  line  between  the  airplane  and 

DME  antenna,  deg 

R  range  from  the  back-azimuth  antenna  to  the  airplane,  m 

rE  radius  of  Earth,  m 

rM  meridional  radius  of  curvature,  m 

rN  normal  radius  of  curvature,  m 

u  unit  vector  perpendicular  to  the  runway  center  line 

vn*ve  north,  east  components  of  navigation  system  velocity  estimate,  knots 

VN,VE  north,  east  components  of  inertial  ground  speed,  knots 

X',Y'  axes  of  orthogonal  coordinate  system  oriented  along  the  runway 
center  line 

•  • 

Xe,Ye  coordinates  of  vector  Ze  transformed  into  the  X',Y'  coordinate 
system 

Ze  vector  of  airplane  estimated  position  from  the  MLS  back-azimuth 

antenna 

Ze  magnitude  of  Ze,  n.  mi. 

ze,N»ze,E  north,  east  components  of  Ze,  n.  mi. 

Z,,,  vector  of  airplane  measured  position  from  the  MLS  back-azimuth 

antenna 

-► 

Zm  magnitude  of  Zm,  n.  mi. 

zm,N'zm,E  north,  east  components  of  Zm,  n.  mi. 

Zr  vector  of  airplane  position  estimated  radially  along  the  measured 

azimuth  angle 

Zr  magnitude  of  position  vector  Zc,  n.  mi. 

zr,N*Zr,E  north,  east  components  of  position  vector  Zr,  n.  mi. 

a  angle  formed  at  the  back-azimuth  antenna  by  the  DME  antenna  and 

the  measured  airplane  position,  deg 


Apn,Ape  north,  east  components  of  position-estimate  error,  n.  mi 
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At 


change  in  time 


AVN,AVE  north,  east  components  of  system  velocity  update,  knots 

A<J>,AX  latitude,  longitude  update  estimates 

n  back-azimuth  angle  relative  to  the  runway  center  line,  deg 

li  relative  angle  between  the  DME  and  the  MLS  back-azimuth  antenna,  deg 

^DME'^DME  latitude,  longitude  of  DME  antenna  location,  deg 

4>e,Ae  latitude,  longitude  of  airplane  position  estimate,  deg 

<|>0,A0  latitude,  longitude  of  back-azimuth  antenna  location,  deg 

ij)r  runway  heading  to  true  north,  deg 

£2  vertical  angular  limit  of  back-azimuth  volumetric  check,  deg 

a)  lateral  angular  limit  of  back-azimuth  volumetric  check,  deg 

Subscripts: 
t  at  time  t 

t-1  at  iteration  time  previous  to  time  t 

Abbreviations: 

DME  distance-measuring  equipment 

IBD  navigation  update  mode:  inertial  velocity,  back  azimuth  and  DME 

IBX  navigation  update  mode:  inertial  velocity,  back  azimuth  alone 

IDD  navigation  update  mode:  inertial  velocity,  dual  EME 

MLS  microwave  landing  system 

NCU  navigation  computer  unit 

PDME  precision  distance-measuring  equipment 

VCR  very-high-frequency  omnidirectional-range  ratio 


AIRPLANE,  EXPERIMENTAL  SYSTEMS,  AND  EQUIPMENT 


Airplane 


The  NASA  test  airplane  is  the  Boeing  737-100  twin-jet  transport  airplane 
shown  in  figure  1 .  The  TCV  airplane  is  used  as  a  research  vehicle  with 
separate  experimental  navigation,  guidance,  flight  control,  and  display 
systems  located  in  a  separate  research  cockpit.  All  of  the  normal  flight 
systems  (flight  control,  navigation,  pressurization,  etc.)  have  been  retained 
in  the  conventional  cockpit  in  a  normal,  functional  state.  This  arrangement 
allows  changes  to  be  made  to  any  of  the  experimental  systems  while  retaining 
the  standard  operational  features  of  the  airplane. 


Experimental  Systems 

The  experimental  research  flight  systems  consist  of  a  triplex  digital 
flight  control  computer  system,  an  electronic  cathode-ray-tube  (CRT)  display 
system,  and  a  digital  navigation  and  guidance  system  (ref.  2)  integrated  into 
a  separate,  two-man-crew,  research  cockpit.  Figure  2  is  a  simplified  func¬ 
tional  block  diagram  of  the  navigation,  guidance,  and  control  process  during 
automatic  path  tracking  (ref.  3) .  The  navigation  and  guidance  calculations 
are  performed  in  a  single  digital  navigation  computer.  Various  navigation 
sensor  signals  (including  inertial-navigation-system  velocities  and  accel¬ 
erations,  true  airspeed,  magnetic  heading,  and  MLS  back -azimuth,  VCJR,  and 
DME  radio  signals)  are  used  in  the  navigation  computer.  This  computer  esti¬ 
mates  the  airplane  position  based  on  combinations  of  the  sensor  inputs. 
Horizontal  path,  vertical  path,  and  thrust  commands  based  on  the  estimated 
position,  velocity,  and  path  tracking  errors  of  the  airplane  are  computed 
and  transferred  to  the  flight  control  computer  system  20  times  per  second. 

The  flight  control  system  then  comnands  the  flight  control-surface  servos. 


Airborne  Data  Acquisition  System 

A  wide-band  magnetic  tape  recorder  onboard  the  airplane  recorded  data  at 
a  rate  of  40  samples  per  second.  These  data  included  93  parameters  describing 
the  airplane  configuration,  attitude,  and  control-surface  activity.  Thirty-two 
additional  channels  were  used  for  recording  parameters,  including  the  latitude 
and  longitude  of  the  navigation  position  estimate  which  was  calculated  by  the 
navigation  computer.  Video  recordings  of  the  electronic  attitude-director  and 
electronic  horizontal-situation  displays  were  also  made  throughout  the  flights. 


Radar  Tracking  Facilities 

Radar  tracking  of  the  airplane  was  provided  by  the  Federal  Aviation 
Administration  Extended  Area  Instrumentation  Radar  (EAIR)  Facility  at  the 
National  Aviation  Facilities  Experimental  Center  (NAFEC)  in  Atlantic  City, 
New  Jersey.  The  tracking  radar  is  a  precision,  C-band,  instrumentation 
radar  system  which  was  operated  in  a  secondary  (beacon- tracking)  mode  during 
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these  flight  tests.  Slant  range,  azimuth  angle,  and  elevation  angle  data 
were  recorded  at  a  10-Hz  sample  rate  on  magnetic  tape.  All  airborne  and 
ground-based  recorded  data  were  time  correlated  for  post-flight  processing 
and  analysis.  Post-flight  processing  of  the  range,  azimuth,  and  elevation 
data  consisted  of  conversion  to  latitude,  longitude,  and  altitude. 

Radar  accuracy,  root  mean  square  (rms) ,  is  stated  (ref.  4)  as  0.15  mrad, 
5.6  m  at  37  040  m  (20  n.  mi.)  in  azimuth  and  elevation  and  6  m  in  range. 


MLS  and  DME 

The  MLS  used  during  these  flight  tests  was  a  test  bed  system  oriented 
for  approaches  to  runway  31  at  NAFEC.  This  system  was  a  time- referenced 
scanning-beam  system  with  a  front-azimuth  beam  width  of  1°,  a  lateral  sweep 
angle  of  ±60°,  and  nominal  vertical  coverage  up  to  20°.  The  glide  slope 
beam  was  1°  wide  and  had  an  azimuthal  coverage  of  approximately  ±60°.  An 
MLS  PDME  with  omnidirectional  coverage  of  approximately  37  040  m  (20  n.  mi.) 
was  used  in  conjunction  with  the  front-azimuth  signals.  The  back-azimuth 
beam  was  3°  wide  and  had  a  lateral  sweep  angle  of  ±40°.  Vertical  coverage 
was  between  1°  and  20°. 

The  Atlantic  City  DME,  normally  used  for  en  route  navigation,  was  used 
during  the  back-azimuth  flight  testing.  This  DME  is  located  at  the  airport 
1671  m  past  the  back-azimuth  antenna  and  229  m  left  of  the  runway  center 
line  (fig.  3). 

When  the  navigation  system  updated  its  position  estimate  with  dual  DME, 
the  Atlantic  City  DME  was  always  used  in  conjunction  with  other  DME  in  the 
area.  The  navigation  system  automatically  selected,  tuned,  and  checked  the 
validity  of  the  other  DME  on  the  basis  of  strength  of  signal  and  station 
location  geometric  characteristics  so  that  good  position  estimates  resulted. 
These  geometric  characteristics  are  described  in  reference  2. 


Switching  Logic  of  Airborne  MLS  Receiver  Antenna 

The  TCV  airplane  was  equipped  with  two  MLS  receiver  antennas  to  ensure 
that  the  greatest  possible  signal  strength  was  delivered  to  the  MLS  receiver 
regardless  of  the  airplane  direction  of  flight,  its  attitude,  or  its  position 
within  the  MLS  signal  coverage.  The  forward  antenna  was  located  on  top  of  the 
fuselage  just  above  the  forward  cockpit.  The  rear  antenna  was  mounted  on 
the  bottom  of  the  fuselage  below  the  airplane  tail.  The  MLS  airborne  receiver 
contained  the  logic  to  determine  which  antenna  was  providing  the  greatest 
signal  strength.  The  receiver  sent  a  pair  of  discrete  logic  signals  to  an 
external  antenna  switch  to  select  the  appropriate  antenna. 

Additional  hardware  logic  in  the  flight  control  computers  forced  selection 
of  the  forward  antenna  while  operating  in  the  MLS  front-course  signal  coverage, 
regardless  of  the  relative  signal  strengths  at  the  forward  and  aft  antennas. 
This  additional  antenna-selection  logic  was  required  because  the  flight  control 
computers  did  not  have  the  necessary  software  for  antenna  switching. 
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EXPERIMENT  DESIGN  AND  PROCEDURES 


The  objective  of  this  flight  experiment  was  to  determine  the  accuracy 
of  the  navigation  position  estimate  when  updated  with  dual  EME  signals,  with 
MLS  back-azimuth  signals,  and  with  a  combination  of  MLS  back-azimuth  and 
DME  signals.  The  flight  and  radar  data  obtained  during  these  tests  were  used 
to  compare  the  accuracies  of  these  update  modes  and  to  determine  the  relative 
errors  due  to  DME  and  back-azimuth  signal  inputs. 

The  position-estimate  error  was  determined  by  calculating  the  vector 
between  the  estimated  position  of  the  airplane  and  its  position  as  tracked 
by  ground-based  radar  (ref.  5) .  The  magnitude  of  this  error  vector  is  the 
position-estimate  error .  The  mean  and  standard  deviation  of  the  navigation 
position-estimate  error  was  calculated  and  used  for  comparison  of  the  three 
update  inodes. 

The  NASA  TCV  airplane  was  flown  along  three,  three-dimensional  flight 
paths  (fig.  3),  each  stored  in  the  navigation  data  base  in  the  navigation 
computer.  Path  1  was  a  5556-m  (3-n.  mi.)  downwind  leg  with  a  180°,  descending 
left  turn  to  a  5556-m  (3-n.  mi.)  final  approach  to  the  runway.  This  path  was 
in  the  MLS  front-course  signal  coverage  and  allowed  the  navigation  computer 
to  update  its  position  estimate  with  the  front-course  signals  just  before 
each  test  run. 

Paths  2  and  3  provided  experience  with  guidance  in  the  area  of  back- 
azimuth  signal  coverage.  These  paths  overlaid  the  last  5556  m  (3  n.  mi.) 
of  path  1  so  that  lateral  guidance  would  be  continuous  during  the  transition 
from  front-course  to  back -course  guidance.  Path  2  extended  18  520  m  (10  n.  mi.) 
from  the  runway  threshold,  straight  out  along  the  runway  center  line.  Path  3 
followed  the  runway  center  line  for  approximately  6482  m  (3.5  n.  mi.),  then 
made  a  45°  left  turn,  and  continued  an  additional  21  298  m  (11.5  n.  mi.). 

Each  test  run  was  started  with  the  pilots  flying  an  instrument  approach 
to  the  runway  via  path  1 .  Instead  of  landing,  a  missed  approach  was  executed, 
and  the  airplane  was  flown  in  the  MLS  back-azimuth  coverage  along  either  path  2 
or  path  3.  upon  reaching  the  end  of  either  of  these  paths,  the  airplane  was 
returned  to  the  beginning  of  path  1  for  another  test  run. 

The  radar  was  calibrated  with  the  airplane  on  a  known  ground  test  point. 
The  radar  tracked  the  airplane  continuously  during  the  flight  from  take-off 
to  landing. 


DESCRIPTION  OF  NAVIGATION  POSITION  ESTIMATE 


Position  Difference 

The  TCV  navigation  computer  was  prog  r  a  timed  to  select  and  tune  two  appro¬ 
priate  DME  and/or  VOR  stations  in  the  vicinity  of  the  airplane.  The  distance 
and/or  azimuth  information  received  from  these  stations  was  used  to  update 
the  previous  position  estimate.  This  position  difference  was  divided  into 
north  (AP(jft)  and  east  (APe,^  components  to  be  used  in  determining  new  posi- 
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tion  and  velocity  estimates.  The  dual  EME  update  mode,  the  primary  position- 
estimate  difference  mode,  was  displayed  as  IDD  on  the  electronic  map  displays 
for  each  of  the  flight  crews. 

Other  sources  of  navigation  information  could  also  be  used  to  determine 
position  differences.  While  operating  in  the  MLS  back-azimuth  signal  cover¬ 
age,  these  position-update  components  were  calculated  with  back-azimuth  and 
DME  information.  This  update  mode  was  shown  on  the  pilots'  electronic  map 
display  as  IBD  (inertial  velocity,  back  azimuth,  DME) .  If  no  DME  was  available, 
the  estimate  mode  was  IBx  (inertial  velocity,  back  azimuth  alone) .  In  subsequent 
sections,  the  equations  used  in  the  IBD  and  IBX  update  modes  are  presented. 

The  equations  used  in  the  IDD  update  mode  are  presented  in  reference  2. 

In  the  TCV  navigation  system  the  same  position-estimate  calculations, 
described  in  the  next  section,  were  used  regardless  of  the  position-update 
mode  used  to  determine  the  north  and  east  components. 

Position  Estimate 

The  first  step  in  the  position-estimate  process  was  to  develop  a  velocity- 
estimate  update  from  the  position  update: 


AvN,t  =  AvN,t-l  +  k2  ApN,t 
AvE,t  =  AvE,t-l  +  k2  ApE,t 


An  updated  velocity  estimate  was  then  made  by  summing  the  velocity  update 
with  ground  speed  obtained  from  the  inertial  navigation  system: 


vN,t  =  AvN,t  +  vN,t 


vE,t  =  avE,t  +  VBft 


A  position  update  in  terms  of  latitude  and  longitude  was  obtained  from 
the  system  velocity  and  the  position-difference  components  as  follows: 


vN,t  At  +  K1  ApN,t 


vE,t  At  +  K1  ApE,t 


This  latitude  and  longitude  position  update  was  based  on  an  oblate-spheroid 
Earth  model  (ref.  2)  by  using  the  following  radii  of  curvature  in  the  meridional 
(north/south)  and  normal  (east/west)  directions: 


rM,  t  *  11  a/p,  t  +  rE(l  -  2F  +  3F  sin2  4>e,t-l> 


rN,t  =[ha/p,t  +  rE(l  +  F  sin2  4>e, t-1  >  1  cos  <t>e,t-l 


The  updated  position  estimate  was  found  by  summing  the  previous  position 
estimate  with  the  position-update  terms: 


*e,t  =  ^e,t-1  +  A<t>t 

^e,t  =  ^e,t-l  +  A^t 

Dual  DME  (IDD)  Position  Estimate 

Dual  DME  updating  is  the  normal  position  update  mode.  The  equations  used 
to  determine  position  difference  while  in  this  mode  are  given  in  reference  2. 


Back-Azimuth/DME  (IBD)  Position  Estimate 

General  solution.-  Figure  4  shows  the  geometry  of  the  MLS  back-azimuth 
antenna,  the  DME  antenna,  the  position  estimate  of  the  airplane,  its  position 
measured  with  back-azimuth  and  DME  information,  and  the  runway  with  an  extended 
center  line.  An  orthogonal  coordinate  system,  with  its  origin  placed  on  the 
back-azimuth  antenna  is  oriented  with  respect  to  true  north.  Since  the  DME 
could  be  selected  either  manually  or  automatically,  the  position-update  equa¬ 
tions  were  derived  so  that  the  DME  was  not  required  to  be  colocated  with  the 
back-azimuth  antenna. 

By  determining  the  relative  geometry  between  the  back-azimuth  antenna,  the 
DME  antenna,  and  the  position  of  the  airplane  measured  with  the  back-azimuth  and 

DME  information,  the  vector  Zmj t  between  the  origin  and  the  measured  position 

of  the  airplane  was  found.  The  vector  between  the  origin  and  the  pre¬ 

viously  estimated  position  of  the  airplane  could  be  found  directly  from  their 
known  latitudes  and  longitudes.  The  position  difference  was  then  found  by  sub- 

-+■  > 
tracting  Zm>t  from 


Calculation  of  the  airplane  position  vector  Z^t*-  Figure  4  shows  that 

Zmjt  is  one  side  of  the  triangle  formed  by  the  origin  (back-azimuth  anlenna), 
the  measured  position  of  the  airplane,  and  the  DME  antenna.  Known  quantities 

used  to  determine  Zm  t  in  this  triangle  include  the  runway  heading  4^ r ,  the 
slant  range  of  the  DME  D^,  the  back-azimuth  angle  relative  to  the  runway  center 
line  n,  the  location  (latitude  and  longitude)  of  the  DME  antenna  'I’dmE'  ^CMF' 
and  the  back-azimuth  antenna  location  4>0,X Q . 

Figure  5  shows  the  angular  geometry  and  distance  A  between  the  origin 
and  the  DME  antenna .  The  length  A  and  the  angle  between  the  DME  antenna 
from  the  north  axis  y  remained  constant  for  a  particular  DME  and  were  cal¬ 
culated  only  once.  If  a  different  DME  was  tuned,  then  A  and  y  were 
recalculated.  The  length  of  A  in  nautical  miles  was  determined  by  vec¬ 
tor  ially  summing  its  components,  a  and  bs 


($o  + 

a  =  <*EME  “  *o>  <6°)  cos  \ - 

b  =  (^eme  "  ^o) (60) 

A  =  /a^  +  b^ 


The  angle  y  was  found  by 


The  angle  ot,  formed  by  side  A  and  vector  Z,,, ^ t  may  vary  continuously 
and  was  calculated  20  times  per  second.  It  ranged  between  0  and  it ,  inclusive, 
and  was  found  from 

at  =  |  2tt  -  y  +  (4>r  -  nt)  | 

If  this  calculated  value  of  at  was  greater  than  it,  then 

at  =  |  2tt  -  [|2tt  -  y  +  (i|»r  -  nt)  |]  | 

The  magnitude  of  the  DME  reading  D^.,  measured  in  the  airplane,  was  the 
slant  range  distance  between  the  ground-based  DME  antenna  and  the  airplane. 
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This  distance  was  slant-range-corrected  to  determine  the  ground  distance  Dt 
between  the  airplane  and  DME 


°t 


sm 


cos 


-1 


Angle  P,  formed  by  the  vector  Zm  t 
and  was  calculated  20  times  per  second . ' 


and  side  D^,  could  vary  continuously 
Prom  the  relation. 


A  °t 

sin  Pt  sin 


angle  P^  was  determined  from 


P^  =  sin-1 


sin  ot+ 


Comparison  of  the  square  of  side  A  with  the  sum  of  the  square  of  side 
and  the  square  of  the  magnitude  of  the  vector  Z6ft_i  (estimated  airplane 
position  vector)  determined  whether  angle  P^  was  obtuse  or  acute.  Hence, 


r 


p*  = 


f  A 

l- 


sin  ■  |  —  sin  a,.) 

\Dt  7 


■< 


/  A 

.  —  1  t  A 

TT  -  Sin  1  I  — 


sin  at  ] 


if  A2  <  Dt2  +  Ze> t-1 2 


if  A2  >  Dt2  +  Ze>t_i 2 


The  magnitude  of  the  estimated  airplane  position  vector  Z6ft_i  was  used  as  an 
approximation  for  Z^.  The  magnitude  of  Ze,t-1  was 


ze,t-l  =  ^J<2e,N,t-1>2  +  (ze,E,t-l>2 
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in  which 


ze,N,t-1  ■  “  4>o)(60) 

/<t>0  +  4>e,t-l\ 

ze,E,t-l  =  ^e,t-l  “  cos  l  ~  J 

-*■ 

The  actual  magnitude  of  Zmf t  was  found  fran 

zm, t  =  A  cos  at  +  Dt  cos  Pt 

-► 

The  north  and  east  components  of  Zmf t  were  calculated  from  the  angle  ipr  -  ht 
between  the  north  axis  and  Zm/t: 

zm,N,t  -  Zm,t  cos  ^r  "  nt> 

ZmjE, t  *  zm, t  sin  <*r  “  Ht) 

Calculation  of  the  position  difference  in  north  and  east  components  APN>t, 

APE>t  —  IBD  update  mode.-  The  position  difference  in  north  and  east  components 

was  now  found  by  subtracting  the  north  component  of  from  the  north 

component  of  Z^  t  and  by  subtracting  the  east  component  of  ze,t-1  from  the 
•¥ 

east  component  of  zmj  t: 

^N, t  =  ze,N,t-1  “  zm,N,t 
^Eft  =  ze,E,t-l  "  zm,E,t 


These  position  difference  components  were  then  used  directly  in  the  navigation 
position-estimate  algorithms. 


Back-Azimuth  Only  (IBX)  Position  Estimate 


General  solution.-  In  the  event  that  a  valid  DME  signal  could  not  be 
obtained,  the  navigation  computer  utilized  the  back-azimuth  signal  and  inertial 
velocity  to  determine  a  position  estimate.  Figure  6  shows  the  geometry  of  the 
back-azimuth  antenna,  the  position  estimate  of  the  airplane,  its  estimated 
position  on  a  measured  back-azimuth  radial,  and  the  runway  with  an  extended 
center  line.  An  orthogonal  coordinate  system,  with  its  origin  at  the  back- 
azimuth  antenna  (<t>0,X0),  is  oriented  in  a  true-north  direction. 

Since  no  EME  information  was  available,  no  radio  position  difference  could 
be  developed  along  a  radial  from  the  back-azimuth  antenna.  However,  inertial 
velocity  was  utilized  in  the  radial  direction  to  supply  inputs  for  a  new  posi¬ 
tion  estimate  in  the  position-estimate  algorithms.  Radio  position  difference 
in  the  IBX  update  mode  was  limited  to  a  direction  perpendicular  to  the  runway 
center  line. 

North  and  east  position  estimates  were  found  in  the  following  manner.  A 

position-difference  vector  DPt  was  found  by  subtracting  the  estimated  position 

-+■  ■> 
vector  from  the  estimated  position  vector  Zr>t  on  the  measured  back- 

azimuth  angle.  The  component  of  DPt  perpendicular  to  the  runway  was  found 
and  broken  into  north  and  east  components.  These  components  were  used  in  the 
navigation  position-estimate  algorithms.  This  process  was  repeated  20  times 
per  second . 

Calculation  of  DPt.-  To  determine  DPt  it  was  necessary  to  calculate  the 


posit ion- estimate  vector  Ze  t_i  in  north  and  east  components  from  the  latitudes 
and  longitudes  of  the  last  position  estimate  and  of  the  back-azimuth  antenna 
locations 


ze,N,t-l  =  <6°)  ^e,t-1  ~  <fr0> 


Ze,E,t-l 


(60)(Xe,t-l 


XQ) 


Ae,  t-1  +  <t>o\ 
cos  y — - — j 


ze,t-l  =  ^e^t-l*1  +  (ze,N,t-1>3 


A  vector  ZCft  of  the  estimated  position  of  the  airplane  along  a  measured 
back-azimuth  radial  was  determined  in  the  following  manner.  Since  no  radio 
updates  could  be  obtained  along  a  radial  from  the  back-azimuth  antenna,  it  was 
assumed  that  the  estimated  radial  distance  from  the  origin  was  correct.  Hence, 
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the  vector  lengths  of  Zr>t  and  Ze,t-1  were  the  same  (directions  could  differ 

— ►  >  “► 

to  obtain  DP^.)  .  The  length  of  ZCjt  and  Ze>t_-j  was 

zr,t  =  ze,  t-1  =  fze,N,t-l)2  +  <ze,E,t-l)2 

->  . 

The  north  and  east  components  of  Zr>t  were  found  from  the  angle  vr  - 
between  the  vector  ZCft  and  the  north  axis: 

zr,N,t  -  zr ,t  cos  Wr  ‘  V 
zrfE,t  =  zr,t  sin  <^r  "  ^ 
zr,t  =  (zr,E,t^  +  ^zr,N,t^3 

■  ► 

DPt  was  found  vec tori ally: 

DPt  =  zr,t  _  ze, t-1 

Calculation  of  the  position  difference  Ap^^.,  ^pE,t  “  IBX  update  mode.- 
The  magnitude  of  the  component  of  DP^-  perpendicular  to  the  runway  center  line 
was  obtained  by  the  dot  product  of  DPt  and  a  unit  vector  u  perpendicular  to 

A 

the  runway  center  line.  The  unit  vector  u,  shown  in  figure  6,  was 
u  =  -cos  (4»r)i  +  sin  (4>r )  j 
The  magnitude  of  DPp,t  was  thus 

DPp,t  -  i*t  •  u 

-  " (zr ,E, t  -  ze,E,t-l>  cos  *r  +  <zr,N,t  '  ze,N,t-1>  sin  $r 

This  provided  the  magnitude  of  the  position  difference.  The  north  and  east 
components  of  the  position  difference  were 
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A*N,t  =  DPp,t  sin  'I'r 


ApE,t  =  “DPp,t  cos  4»r 

These  position-difference  components  were  used  directly  in  the  navigation 
position-estimate  algorithms. 


Back -Azimuth/Airplane  Position  Validity  Check 

Since  random  DME  selection,  automatic  frequency  tuning,  and  other  means 
of  automatic  software  control  were  utilized  in  the  navigation  computer,  a  check 
had  to  be  made  to  ensure  that  the  navigation  computer  was  using  the  appropriate 
navigation  data  for  the  area  in  which  the  airplane  was  flying.  Obviously,  if 
improper  navigation  data  were  being  utilized,  position  estimates  based  on  those 
data  would  be  inaccurate  and  should  not  be  used. 

To  preclude  the  use  of  the  wrong  MLS  back-azimuth  navigation  data,  a  volu¬ 
metric  geometric  check  was  made  to  determine  whether  the  position  estimate  of 
the  airplane  was  within  the  MLS  back -azimuth  boundaries.  If  the  airplane  was 
not  within  these  boundaries,  then  MLS  updating  was  inhibited.  These  boundaries 
included  lateral  and  radial  limits  of  coverage  as  shown  in  figure  7.  The 
lateral  limit  of  the  back-azimuth  angle  of  coverage  is  ±to  from  the  runway 
center  line.  Radial  limits  required  that  the  airplane  be  within  a  distance  L 
in  meters  (n.  mi)  of  the  back-azimuth  antenna,  but  not  closer  than  a  distance  M 
in  meters  (n.  mi) .  A  vertical  angle  of  coverage  limit  ft  was  measured  from 
the  back -azimuth  antenna. 

A  new  orthogonal  coordinate  system,  with  its  origin  located  at  the  back- 
azimuth  antenna  and  its  X'-axis  parallel  to  the  runway  center  line,  was  used 
to  make  the  geometric  check.  The  north  and  east  components  of  the  estimated 

position  vector  Ze,t-1  were  transformed  into  the  new  X',Y'  coordinate  system 
as  follows: 


xe,t  "  ze,N, t-1  cos  ’I’r  +  ze,E,t-l  sin  ^r 
Ye,t  =  ze,N,t-l  sin  *r  “  ze,E,t-l  cos  <Pr 
The  lateral  azimuth  check  was  then 


-Xgft  tan  o)  S  Ye,t  “  xe,t  tan  u 
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The  radial  check  was 


M  S 


COS  ht 


£  L 


The  vertical  check  was 


ha/p,t  =  ho  +  xe,t  tan  ^ 


For  these  tests,  the  MLS  back-azimuth  coverage  limits  were 


Lateral .  ^  *  ±42° 

Vertical .  ft  =  20° 

Inner  radial .  M  =  3048  m 

Outer  radial .  L  *  37  040  m  (20  n.  mi.) 


RESULTS  AND  DISCUSSION 

Flight  data  were  gathered  on  a  total  of  43  test  runs  along  path  2  and 
path  3  and  on  the  return  flight  to  the  start  of  path  1 .  Position-estimate 
updates  were  accomplished  on  27  update  segments  in  the  IBD  mode,  7  update  seg¬ 
ments  in  the  IBX  mode,  and  30  update  segments  in  the  IDD  mode.  Typically,  seg¬ 
ments  lasted  between  1  and  4  minutes.  Position  estimates  from  the  navigation 
system  in  each  update  mode  were  compared  with  radar  data  to  obtain  the  position- 
estimate  error  vector.  The  mean  and  standard  deviation  of  the  magnitude  and 
direction  of  this  vector  were  determined  for  each  update  mode. 


IBD  Update  Mode 


In  the  IBD  mode,  the  magnitude  of  the  position-estimate  error  vector  had 
a  mean  of  124.3  m  and  a  standard  deviation  of  60  m.  The  mean  direction  of  the 
error  vector  was  298.1°  with  a  standard  deviation  of  35.7°. 

Figure  8  shows  how  the  position-estimate  error  vector  is  divided  into  two 
components,  so  that  the  errors  due  to  the  back-azimuth  signal  input  and  by  the 
DME  signal  input  can  be  separated.  The  first  component  ep  is  perpendicular 
to  the  back-azimuth  radial  on  which  the  airplane  is  located.  The  ep  component 
is  influenced  by  the  back-azimuth  signal  input.  The  second  component  e^  is 
a  radial  error  along  the  back-azimuth  radial  on  which  the  airplane  is  located. 

The  eg  component  is  influenced  primarily  by  the  DME  signal  input  and  to  a  lesser 
extent  by  the  back-azimuth  signal.  The  e^  component  increases  linearly  as 
a  function  of  the  range  R  from  the  back-azimuth  antenna  to  the  airplane  if 
the  azimuth  bearing  error  eg  (fig.  8)  is  constant  (i.e.,  ep  =  Reg). 
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Figure  9  shows  the  ep  component  as  a  function  of  R  for  selected  points 
from  flights  along  the  zero  back -azimuth  radial.  These  data  are  grouped  about 
a  straight  line  with  a  slope  representing  a  constant  eg  =  0.05°.  To  determine 
how  the  ep  component  varied  on  other  back-azimuth  radials,  data  were  used  from 
flights  along  two  other  path  segments.  The  first  path  segment  was  the  last  leg 
of  path  3  which  diverged  from  the  runway  center  line  at  a  45°  angle.  The  second 
path  segment  was  parallel  to  the  runway  center  line,  but  offset  by  6100  m.  This 
path  segment  was  flown  back  toward  the  runway.  These  paths  provided  data  at 
various  ranges  for  all  back-azimuth  radials  on  the  left-hand  side  of  the  azimuth 
coverage.  (Flight  during  the  test  period  was  confined  to  this  side  of  MLS  cov¬ 
erage  due  to  air-traff ic-contiol  considerations.) 

Figure  10  shows  the  ep  component  as  a  function  of  range  from  the  back- 
azimuth  antenna  during  flight  along  these  path  segments.  The  data  lie  above 
the  line  representing  a  constant  eg  =  0.05°.  Therefore,  eg  was  larger  for 
azimuth  radials  which  were  not  on  the  runway  center  line. 

To  determine  how  position  error  varied  for  different  back-azimuth  radials, 
the  ep  component  from  all  flights  (including  those  along  the  zero-azimuth 
radial)  was  normalized  for  range  effects  and  plotted  in  the  form  of  eg  as  a 
function  of  the  back-azimuth  radial.  Figure  11  shows  that  these  data  were 
grouped  about  a  straight  line  drawn  from  eg  =  0.05°  at  the  zero-azimuth 
radial  to  eg  =  0.28°  at  the  40°  radial.  The  straight  line  in  figure  11  shows 
the  expected  eg  error  as  a  function  of  the  back-azimuth  radial.  An  expected 
ep  component  can  be  found  by  multiplying  the  range  from  the  back-azimuth 
antenna  and  the  appropriate  eg  obtained  from  the  straight  line. 

The  radial  component  of  error  eR  was  plotted  as  a  function  of  range 
from  the  back-azimuth  antenna  in  figure  12.  The  eR  components  for  flights 
on  both  path  2  and  path  3  are  randomly  grouped  about  a  mean  value  of  120  m. 

The  figure  shows  that  the  eR  component  did  not  vary  significantly  with  range. 

The  ei*  component  contributed  the  most  to  the  navigation  position- 
estimate  error.  The  ep  component  contributed  a  smaller  error  that  varied 
as  a  function  of  the  back-azimuth  radial  and  the  range  from  the  back-azimuth 
antenna.  The  expected  ep  component  was  shown  to  be  less  than  32  m  within 
a  range  of  37  040  m  (20  n.  mi.)  from  the  back-azimuth  antenna.  For  operations 
off  the  zero-azimuth  radial,  the  ep  component  will  become  as  large  as  the 
eR  component  only  when  the  airplane  is  close  to  the  lateral  limits  of  the 
back-azimuth  coverage  and  at  large  distances  from  the  back-azimuth  antenna. 

The  DME  ground  station  used  during  these  flight  tests  was  a  standard 
DME  station  used  for  normal  en  route  navigation.  The  expected  accuracy  of 
these  stations  is  about  185  m  (ref.  6).  A  reduction  of  the  DME  ground- 
station  error,  which  could  be  obtained  with  a  PDME,  could  reduce  the  radial 
component  of  error  and  result  in  an  improved  navigation  position  estimate  dur¬ 
ing  the  1BD  update  mode. 
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IBX  Update  Mode 


In  the  IBX  update  mode,  the  magnitude  of  the  position-estimate  error 
vector  had  a  mean  of  145.3  m  and  a  standard  deviation  of  71.3  m.  The  mean 
and  standard  deviation  for  the  direction  of  the  position-estimate  error 
vector  was  not  summarized  since  radio  position  updating  does  not  occur 
omnidirectionally  in  this  mode. 

Both  the  magnitude  and  direction  of  the  position-estimate  error  vector 
were  strongly  dependent  upon  the  initial  error  when  the  IBX  update  mode 
started.  Corrections  to  the  position  estimate  were  made  only  in  a  direction 
perpendicular  to  the  back-azimuth  radial  on  which  the  airplane  was  located. 

No  corrections  could  be  made  in  a  radial  direction  since  DME  signal  inputs 
were  not  available. 

Figure  13  shows  the  ep  and  e^  error  components  for  a  flight  along  the 
zero-azimuth  radial.  The  IBX  update  mode  was  started  with  initial  navigation 
position-estimate  errors  of  ep  =  140  m  and  e^  »  118  m.  The  ep  component 
was  reduced  because  of  back -azimuth  signal  inputs  but  the  eR  component 
remained  virtually  unchanged  since  no  DME  signal  inputs  were  available. 

Figure  14  shows  the  ep  component  as  a  function  of  range  when  the  air¬ 
plane  was  flown  along  path  5  and  back  to  the  start  of  path  1 .  The  results 
were  similar  to  those  obtained  in  the  IBD  update  modes  eQ  =  0.05°  on  the 
zero-azimuth  bearing  and  efl  £  0.25°  on  the  return  leg. 


IDD  Update  Mode 

In  the  IDD  update  mode,  the  position-estimate  error  vector  was  random 
in  direction  and  had  a  mean  magnitude  of  183.5  m  and  a  standard  deviation  of 
81.1  m.  These  results  are  within  the  185-m  accuracy  tolerance  for  DME  ground 
stations  and  are  similar  to  past  experience  with  the  dual  DME  update  mode 
(ref.  5) .  It  is  expected  that  the  use  of  PDME  inputs  also  could  improve  the 
position-estimate  accuracy  in  the  IDD  update  mode. 


CONCLUDING  REMARKS 

Flight  tests  with  the  NASA  Terminal  Configured  Vehicle  (TCV)  Boeing 
737-100  airplane  in  the  time- referenced  scanning-beam  (TRSB)  microwave 
landing-system  (MLS)  environment  have  shown  that  the  accuracy  of  position  esti¬ 
mates  based  on  dual  distance-measur ing-equipment  (DME)  inputs  can  be  improved 
by  using  MLS  back-azimuth  radio  signals.  Position  estimates  from  a  combination 
of  back-azimuth  and  DME  signals  achieved  the  best  position-estimate  accuracy. 
The  next  best  accuracy  was  obtained  from  back-azimuth  signals  alone.  The  least 
accurate  update  mode  of  the  three  tested  was  the  dual  DME  mode. 

The  position-estimate  accuracy  was  established  by  computing  an  error 
vector  between  the  position  of  the  airplane  defined  by  ground-based  radar 
tracking  and  the  airplane  position  estimated  by  the  airplane  navigation  system. 
This  error  vector  was  divided  into  two  components:  a  component  perpendicular 


to  the  back-azimuth  bearing,  influenced  by  back-azimuth  signal  inputs;  and  a 
component  parallel  to  the  back-azimuth  bearings,  influenced  primarily  by  the 
DME  signal  inputs. 

For  updates  using  both  DME  and  back-azimuth  signal  inputs,  the  radial 
component  of  error  contributed  the  largest  portion  to  the  error  vector.  The 
radial  component  of  error  was  generally  constant  as  a  function  of  range  from 
the  back-azimuth  antenna  and  was  within  accuracy  tolerances  for  this  DME 
station.  The  DME  ground  station  used  during  these  tests  was  a  standard 
en  route  navigation  station.  The  use  of  a  precision  DME  ground  station 
probably  could  improve  the  accuracy  of  the  position  estimate. 

The  perpendicular  component  of  error  was  influenced  both  by  range  from 
the  back-azimuth  antenna  and  by  the  particular  back-azimuth  radial  on  which 
the  airplane  was  located. 

For  updates  which  used  only  the  back-azimuth  signal  inputs,  position- 
estimate  error  was  strongly  dependent  upon  the  initial  position-estimate 
error  since  corrections  were  made  only  in  the  perpendicular  direction.  The 
perpendicular  component  of  error  showed  accuracy  characteristics  similar  to 
those  observed  during  updates  using  both  DME  and  back-azimuth  inputs. 

For  updates  which  used  dual  DME,  position-estimate  error  was  within  the 
accuracy  tolerance  specified  for  en  route  DME  and  was  similar  to  previous 
experience  with  dual  DME  updating.  The  accuracy  of  the  position  estimate 
could  be  improved  with  the  use  of  precision  DME  signal  inputs. 
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Figure  1 NASA  terminal  configured  vehicle  (TCV)  research  airplane. 


Flight  test  paths  within  the  MLS  front-  and  back-azimuth  signal  coverage 
Headings  are  relative  to  true  north. 


position  estimate, 


'igure  7.-  Lateral  and  radial  MLS  back-azimuth  geometry  limits:  U),  l,  and 


eg  are  positive  as  shown. 


Perpendicular  component  of  position  estimate  error  on  various  back-azimuth 
radials  up  to  40°;  IBO  update  mode. 


Perpendicular  component  of  position-estimate  error  normalized  for  range 
a  function  of  the  back-azimuth  radial;  IBD  update  mode. 


Figure  13.-  The  ep  and  eR  components  of  position-estimate  error  from  an 

initial  error;  1BX  update  mode. 
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